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[i] Nongyrotropic plasma distribution functions can be formed in regions of space where guiding 
center motion breaks down as a result of strongly curved and weak ambient magnetic fields. Such are 
the conditions near the current sheet in the Earth’s middle and distant magnetotail, where observations 
of nongyrotropic ion distributions have been made {Frank et ai, 1994]. Here a systematic parameter 
study of nongyrotropic proton distributions using electromagnetic hybrid simulations is made. We 
model the observed nongyrotropic distributions by removing a number of arc length segments from a 
cold ring distribution and find significant differences with the results of simulations that initially have 
a gyrotropic ring distribution. Model nongyrotropic distributions with initially small perpendicular 
thermalization produce growing fluctuations that diffuse the ions into a stable Maxwell ian-1 ike 
distribution within a few proton gyro periods. The growing waves produced by nongyrotropic 
distributions are similar to the electromagnetic proton cyclotron waves produced by a gyrotropic 
proton ring distribution in that they propagate parallel to the background magnetic field and occur at 
frequencies on the order of the proton gyrofrequency. The maximum energy of the fluctuating 
magnetic field increases as the initial proton distribution is made more nongyrotropic, that is, more 
highly bunched in perpendicular velocity space. This increase can be as much as twice the energy 
produced in the gyrotropic case. INDE^' Terms: 7867 Space Plasma Physics: Wave/particle 
interactions, 7843 Space Plasma Physics: Numerical simulation studies, 2772 Magnetospheric 
Physics: Plasma weaves and instabilities, 2764 Magnetospheric Physics: Plasma sheet; KEYWORDS: 
Plasma instabilities, nongyrotropic distributions, hybrid simulations 


1. Introduction 

[2] The vast majority of studies concerning linear and nonlinear 
plasma instabilities in the magnetosphere involve particle distribu- 
tions that are gyrotropic. Part of the reason for this is that in many 
regions of the magnetosphere, g>TOtropic, Maxwellian-like distri- 
butions are commonly observed and pro\ ide an adequate descrip- 
tion of the plasma. Recent observations and theory have led to an 
interest in understanding non-Maxwellian particle distribution 
functions that can be classified as nongyrotropic. A nongyrotropic 
plasma is characterized by a velocity distribution that satisfies 


where f is the three-dimensional distribution function in cylindrical 
velocity coordinates, the velocity coordinates are given by Vj = v>, 
= \\^ + v-‘, and C> = tan“'(iyv,.) is the gyrophase angle in the 
perpendicular velocity plane. The directions parallel and perpendi- 
cular are with respect to the background magnetic field (indicated 
by Bo). 

[3] Recent observations by the Galileo and Geotail satellites in 
the Earth’s magnetotail have established that ion distribution func- 
tions can have non-Maxwellian, nongyrotropic features [Frank et 
al . , 1 994; Saito et al . , 1 994). In fact, the Galileo data led Frank et al. 
[ 1 994] to conclude that distributions in the intermediate {x ^ 35 Re) 
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and distant magnetotail caimot accurately be described in terms of 
convecting, quasi -Maxwellian distributions. 

[ 4 ] The Geotail spacecraft also observed ring-like proton 
velocity distributions that are somewhat nongyrotropic in the 
region between the plasma sheet and the l^be in the distant tail 
[Saito et al., 1994]. Cold proton beams observed in the lobe are 
believed to be the source of the nongyrotropic protons. These 
beams can be accelerated by thennoelectric fields associated 
with slow mode shocks. If the spatial scale of the electric field 
is smaller than the ion gyro radius, the ions can form a ring- 
shaped distribution [Sai/o et al., 1994]. The perpendicular (to 
Bo) speeds are observed to be on the order of a few' hundred 
kilometers per second (for reference, the AlfVen speed is '^23(X) 
km/s) and ring proton to total proton density' ratios ranging from 
0.20 to 1.00. During a 1 -month period the Geotail spacecraft 
spent 80 hours in the area of the plasma sheet-lobe boundary. 
During that time, 25 examples of ring-like distributions with 
varying degrees of nongyrotropicity were observed. Figure Ic of 
Saito et al. [1994] shows such a distribution. 

[ 5 ] Large-scale modeling studies have also supported the idea 
that naturally occurring nongyrotropic-like distributions may be an 
in^x>rtant feature in the magnetosphere. For example, using large- 
scale kinetic (LSK) calculations, Ashour-Abdalla et al. [1993a] 
showed that azimuthally nonsymmctric distributions can form in 
the central plasma sheet (CPS) that can increase the free energy to 
levels above those of the standard loss cone distributions. LSK 
results have also shown that Maxwellian ion distributions from the 
mantle can lead to the formation of a thin current sheet in which the 
ion pressure tensor has significant off-diagonal terms, indicative of 
the presence of nongyrotropic distributions [Ashour-Abdalla et al., 
1994], These nongyrotropic distributions formed by ion nonadia- 
batic motion in the magnetotail are in good qualitative agreement 
with the main features of the Galileo observations [Ashour-Abdalla 
et al., 1996]. Using a centrifugal impulse model, Delcourt and 
Martin [1994, 1999] also demonstrated that nonadiabatic particle 
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morion m the magnetotail could lead to gyrophase bunchmg when 
the particles do not consen'e their magnenc moment near the 
™g^etotail midplane, where local magnetic field ^ ^ 

the order of the Larmor radius. Particle trajectory studies have » “ 
shown that off-diagonal terms in the pressure tensor of the 
magnetotail are essential for maintaining the pressure bailee 
and the tail structure of a thin current sheet at distances >20 Rb 
from Earth [Ashour-Ahdalla el ai. 1993b], These terms are a d.r^ 
consequence of nongyrotropic distributions and are nonexistent fo 
a S^nopic plasmlBecause of their lack of self-consistency, 
hoover, test particle and LSK methods cannot be used to inves- 
tigate the effects of wave particle interactions or other collective 
effects that could be caused by nongyrotropic distributions. 

[6] The question that arises is whether the nongyrott^ic nature 
of these distributions represents a significant source of free enCT^ 
that can dnve plasma instabilities. If so. this >mply 

nongyrotropic aspects of a velocity distnbution should he ‘n^lud^ 
in th^retical and simulation models applied 

other regions of space where nongyrotropic distnbutions may be 

The relevance of studies of nongyrotropic ion distributions 
is not limited to the magnetotail applications. Observations of 
nongyrotropic ions have been made by the Giotto spacecraft ne 
the wmet P/Grigg-Skjellerup (P'G-S), where 
largest near the comet, have been observed up to several particles 
per cubic centimeter [Coarcs cr a/., 1993]. Perpendicular spet^^e 
on the order of 400 km/s. In a study of low-frequency el^o- 
magnetic waves observed near P/G-S, Neubauer et { ] 

found waves with a peculiar character near to, but upstream of. 
the cometary bow shock and suggested that the unusual ch^cter- 
istics may be due to either the nongyrotropic ions swn there or 
nonlinear^effects. These waves are left-hand polmz^, propagate 
along the background magnetic field (within ±10 ), have frequen^ 
cies close to the local ion cyclotron frequency, and are observed m 
conjunction with the nongyrotropic ions. The unique asp^t of 
these waves is that they are nearly sinusoidal with sharp pe^ and 
loaded troughs Neubauer et ai, 1993 Figuros 10a and 

10b] They note, however, that “much more theoretical work k 
needed to assess the importance 

wave fields” (p. 20,950). Mote recently, Astudilh et a/. [1990] 
used data from the PHOBOS II mission to Mars to provide clear 
evidence of nongyrotropic alpha particles and protons in the solar 

wind a.t ~1.3 AU. . 

[si The existence of nongyrotropic ions in space plasmas is 

established. While the quasi-linear and nonlinear 
these distnbutions have only recently been explored in the litCT- 
ature, studies concerning the linear stability of 
distnbutions have been conducted for several years. Suthr, [1965] 
was the first to derive the dispersion relation for parallel propagat 
mg. electromagnetic instabilities dnven by a nongyrotropic ptema^ 
ffedericks [1975] performed a Imear stability analysis 
waves driven by protons that had been partially gyrophase bunched 
by being reflected off the Earth's bow shock. In subsequent ye^, 
sLlarLdies were conducted on the linear dispersion of parallel 
propagating electromagnetic waves driven by 
ffinncfl et ai, 1992. 1993a] and electrons [Fruchtman and Fned- 
land 1983; Freund et ai, 1987]. Linear theory was appli^ by 
Mo^chmann and Classmeier [1993] to investigate the ‘nstabilmw 
generated by nongyrotropic cometary ions, and Cao et ai [19 , 

1998] studied the linear properties of model nongyrotropic ion 
distributions found near comet P/Grigg-Skjellerup. 

[9] Bnnca and Romeiras [1998] did a 
statiLarv nongyrotropic distnbutions. that is, those distribution 
Scan be mamtained by particle sources and sinks winch can 
used to model the plasma of cometary environments. 
differs from ours in that we will smdy simulations of homoge- 
neous time-varying nongyrotropic distnbutions wi* no sources or 
sinks Brinca and Romeiras [1998] showed that adding a stationary 


nongyrotropic distribution to an otherwise stable p asma po^k- 
tion^Sld bring about linear instability m parallel eigenmod _ 
Brinca [1999] smdied the linear stability of both station^ ^ 
tiSviJying iongyrotropic distributions and showed that th^ 

OT can dcsrabilize a plasma and enhance the grov^ of 
instaWities driven by gyrotropic distnbutions. The 
of both stationary and time-varying 

was also studied by Motschmann and Glassmeter ['998]^d 
Motschmann et ai [1999] for a particular parameter regime^ "^ey 
showed that stationary nongyrotropic distnbutions dn 
perpendicular and parallel propagating waves unstable, whi e tin« 
^rionary nongyrotropic distributions led pnmanly to an instability 
oTTrStrurwaves. Burtnskaya et ai [1994] investigated 
stab^S^properties of electrostatic waves in the presence of ion 
distributions with velocity space holes. Thwe holes ^ “ 
way of modeling the nongyrotropic ions produced by chaotic ion 

acceleration in the neutral sheet. » 

fiol Only a few studies have used simulations to investigate the 

nonlinear evolution of nongyrotropic ions. 

[1993] found that shock-produced gyrophase bimched tons gen- 
erated parallel propagating, low-frequency 

high-fit^uency whistler waves in a one-dimensional hybnd sxmu- 

larion. 2 . some cases the low-frequency waves 

gyrophase bunching through particle trapping, and in others they 

scattered the distnbution toward gyrotropy. 

fill Brinca et ai [1993b] used a one-dimensional electromag- 
netic hybrid code to compare instability growth dnven by aniso- 
rtopic gyrotropic protons versus growth by anisotropic, 
nongy^wpic protons. They made the protons severely nong^ 
trophy restricting their initial gyrophase angles to <t> - 
gro^ of the fluctuating magnetic field energy density was 
S2)ccur faster for the nongyrotropic case; howev^ the 
level of the magnetic energy was about the same. The investigators 
interpreted this to mean that both the gyrotropic and nongyrotropic 
cases were driven by the same free energy, that is, the temperat^ 
anisotropy. Another difference between the two ^ 

^gnetifenergy, which tended to oscillate in the nongY^troP^ 
case as energy was exchanged between the 

protons and the electromagnetic fields. Last, they found that tihe 
Nongyrotropic nature of the protons persisted until several gyro 

’^[I 2 ?^£cTm^!ITr a/. [1997] used a two-dimensional hybrid 
code to study nongyrotropic ions in association with comets by 
injecting spatially homogeneous alpha particles 
thi background magnetic field for a time shorter than m ion 
gyro period. This process produced nongyrotropic distnbuti^ 
*at in velocity space look like rings with two arc length 
removed. The main findings concerned two ^ 

diffusion. Radial diffusion occurred on timescales 
of 10 Up-', where Q, is the proton gyrofrequency. Diffus on m 
the azimuthal direction depended on the size of the ^ 
removed from the ring distribution and could be “ 
diffusion There was no discussion of saturation proces^ 
Recently Cao et ai [2000] used a one-dimensional 
SSfc hyrrid CO* i s Jy *c pickop P-OCCS. and d.lta* 
of Mwbom ions in solar wind plasmas. The p^Uel 
electromagnetic ion cyclotron ring instability dominated and W 
to both phase and pitch angle diffusion on a timescale 
order of a few proton gyro penods. Phase angle scattenng of the 
ions occurred slightly faster than pitch angle scattenng and 

much faster than energy diffusion. -lertm- 

[ 13 ] Brinca et ai [1998] considered one-dimensional elw^^ 
static and electromagnetic particle in cell simulations o 
varying nongyrotropic distribution functions. They found that 
el^trostatic cyclotron harmonic waves could be dnven unstabl 
as ^electromagnetic ordinary mode. The perpend.cuU 

instabilities resulted in ^phase diffusion that accompanies the 
evolution toward saturation. 
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fl 4 l The magnetospheric observations highlight Ae ® 

thorough study of nongyrotropic distributions and their 
consequences in order to deteimine their role in magnetospheric 
Dhvsics The importance of nongyrotropic distributions as a source 
of^ energy to drive instabilities is still not clear, nor is it cl^ i 
fte^m,?ic plasma particle picture of 

man^eoretical studies needs to be modified to include these 
effects This paper will address these issues by cairying a 
sysTe^atic paL^tne simulation study of nongyrotropic distnbu- 
tion functions The goal is to contribute to the understanding of the 

d»,ons »d hdp fo™ . 

complete picture of the magnetotail. 

[i 5 l In section 2 the hybrid simulation code us^ for all of the 
simulations in this paper is described We also discuss in more 
detail the observational foundations for the 
distribution used in this study. Section 3 includes a detail^ 
^alysis of a simulation using a ring distnbution and chares 
these results to several other simulations using different types of 

velt^ space diffusion, and individual electric field and power 
ISI charactenstics are discussed. This section also present 
r^ults of the growing fluemating field dependence on the ini la 
values of the ^rpendicular beam speed, the nongyrotropic proton 
distribution density, and the shape of the nong^otropic distnbution 
in velocity space. A detailed discussion of the saturatiori mecha 
nism involved in the growing waves is also ^ 

section Finally, properties of model nongyrotropic distributions 
that may cause the plasma to maintain stability even though it is 
highlv L-Maxwellian are presented Section 4 summanzes om 
coLl'usions and discusses the relevance of these results to local 
and global magnetospheric physics. 


1 Inless Otherwise Noted in the Text 

Parameter 

Value 

System length 

Grid size 

Time step 

Total run time 

Total number of particles 

Nongyrotropic plasma densiw 

Perpendicular velocity 

Electron beta (3* 

Initial proton beta 3],p 

L, = Ly = 32 c/jp, 

Ax = Ay = c/uop, 

Af = O.OSfip' 

25 < IV < 100 
819.200 < V < 1.310,720 

0.01 < Elng Hc < 10 

1.0 < Vj.p/VA < 5,0 

1.0 

lO-*^ 3,,p< 10 


2. Hybrid Code and Model Distribution 

[ 16 ] The simulations presented in this smdy are c^ed out 
Jng a two-dimensional electromagnetic hybnd code. This code 
w^'onginally written by Dan Winske and is^ed on lu one^ 
dimensional code described by fVinske and Omidi ] 

mam advantage of a hybnd code is the great economy of comput- 
ine time achieved by treating the electrons as a inassless fluid. A 
fiiflv kinetic description of protons is used, so the foil ^ge of ion 
Satial and temporal scales is resolved. The code also sets foe 
d^snlacement current equal to zero in Maxwell s equations. Such a 
scheme is valid for foe smdy of low-frequency e^ctroma^etic 
orocesses such as those that we are interested in. The simulation 
box is periodic in the two spatial dimensions (x-v)^each having 
i rrtVi / = / = where c is the speed of light and 

foe total proton plasma frequency^The system 
Tenefo wi^ chosen so that it adequately resolves foe dominant 
/no mfvie*; as determined by the linear theory' dispersion 
SfotiOTfor a proton ring distribution [Convery 1997] 

[ 17 ] Fluctuating fields are allowed to propagate fi^ly m the (x- 
V) plane. There is a uniform ambient magnetic field in foe x 
dirLon The simulation is earned out in foe reference ffome 
there foe electrons are at rest. Complete simulation parameters are 

protons are loaded uniformly in space. 
grourld protons (if any) have density T 
dunensional. non-drifting Maxwellian velocity 
total proton density' /ip satisfies np = - «b.ck - "c- 

foe nongyrotropic proton density and n. is the total electron 
density ^e nongyrotropic proton distnbution is ^ 

the parallel direction and is modeled by a cold nng i ^ 
perpendicular (vv - v^) velocity plane. The nng distnbution is 
STmodified b^ removing one or more 

various sizes to create a nongyrotropic distnbution. The tou 
TmuLT wfofo of these sections determines foe nongyrotropic 


nature of the distribution at any given time. In section 3. a 
parameter <I> is defined that quantifies the non^otro^c ^pec 
for any given distribution. The parallel Protoi^'as™ betaj^ _ 

Stw T«JB ^ = 10■^ unless otherwise rioted. The Alfoen speed 

■■ I'ijr sf » 

closely resembles foe observations of nongyrotropic particles in fo 
eS's mS^etotail [e.g.. Sdim et ai. 1994]. We also benefit from 
foe added advantage that the simulation results using this mod 
cTn compared to results from simulations using the well- 
understood ifog distribution. This proves to be very helpfol in 
interpreting the simulation results and in isolating uniquely no 
gyi^pic^effects. Both to check our work and to 
S the physics, we contrast our nongyrotropic results with foe 

ring case results throughout this paper, „.;K,.tir.n 

[ 20 ] Although the linear theory for a nongyrotropic 
is not derived here, it is worth discussing the validity of a model 
nonSwVc distribution in terms of linear Vl^v foeoo' m 
antid^tion of future linear theory ’Rations of no^trop^ 

distributions. A nongyrotropic distnbution fun^on ^ 

arido 4 0 where © is foe azimuthal angle in foe velocity plane 
pei^Sidlcular to the ambient magnetic field. Therefore m^der ro 
satisfy the Vlasov equation with no source or smk terms, 
equilforium nongyrotropic distribution must also be e^r time 
dLendent {dfidt / 0) or spatially inhomogMeous ^ ° 

Sn 1965], where i is the three-dimensional position vwtor^h 
is important to note that in a plasma environment wifo local 

particfr sources and sinks, an equilibn^ nongyrotropic fostnbu^ 

iion that is both time stationary and homogeneous be a 
solution to foe Vlasov equation in an open phase space. Alfooug 
noiM here, such a model is appropriate when foe pl^ma 
population includes pickup ions, such as in cometary 
ments (see references in section 1). For foe ptema ^vironrnent 
snidied here (no particle sources or sinks), this implies that 
nongyrotropic distribution must have a functional «‘®^"‘^®"^® 
f ” = /■ ” (V, V . o - ill) for the time-dependent case 

or for the spatialW dep^dlmt case/„/ = /„g‘’ ('„• 'V. o - 
Our model, bunched in velocity space and spatially ^ 

example of the time-dependent form and is therefore a 
solution Although not discussed here, an example of nongyro- 
tropic particle distributions showing spatial bunching was given 
bv'fos^n el ai [1981], who observed ion distributions in foe 
for«hock region displaying gyrophase organization with a spatial 

srale less than an ion gyro radius. 

[■’ 1 ] Motschmann el ai [1997] have done a simulation smdy 
with a two-dimensional hybnd code using a nongyrotropic model 
somewhat similar to what is used here The pnmaiy differences 
with this work are that their work focused on 

cometary environments, specifically, comet P/Gngg-Skjellerup 
wh^e nongyrotropic ions are alpha particles ^d foe ratio of 
S^sity to^l ion density' is of the order of a few percent We 
assoc, .itb «.= Ea,* - ; 

whL nongyrotropic proton densities can be as much as 100/« 
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3. Results of Hybrid Simulations 

3 1 Definition of Parameters 

;:xS 

,e™. of .bo foiling fou, 

parameters: 

Parallel proton beta 


$(,)s [ 

J 


n„g at t = 0 satisfies 


1-K 


0)d<^- 


U„ U«„„ wuh The of 

throughout this paper ^ 1. The top 

distnbutions with diffenmi are shown m 

™dLTanl.'V S 

r rS^inTaS K = .V Staf 

bunches JJ'fnto'LcLt bunch number or spt^ 

rr — > ^ 

total azimuthal two example distributioo! 

simulation results produced by *«e two « P ^ 

imply that should be modified m ^ e« 

other characteristics (bunch "umher an ^ ^ 

distinguish nongyrotropic distn ^ ‘ ^ j 24 The bottor 

.hov^s a two-bunch distribution with 
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Figure 2. 

proton gyrotropic nng 
direction and by fip/ = 
instability. 


Velocitv space time senes in the plane perpendicular to the background magnetic field for an initially cold 
itroDic nnE. Here i ^„'va = 10, and n, = 1-0. The nng distribution rapidly diffus^ m *e radial 
3. the distribution has lost most of the positive gradient [df d\\_ > 0) which drives the 


middle and bottom right panels both shov distributions that 
have 4>o = L5; however, the last distribution (bottom nght 
panel) has a significant thermal spread in the perpendicular 
velocity. By comparing these last two cases, we seek to 
determine the effects of a nongyrotropic distribution with a 
perpendicular thermal spread. The thermal distribution is made 
such that it has no dfd\ ^ > 0 slope, so if it were a gyrotropic 
distribution, it would be stable. We will determine whether or 
not making this thermal distribution nongyrotropic (i.e.. such 


that the distribution has dfdo > 0) by removing two large arc 
length segments will cause it to be unstable. 

3,2. Individual Runs 

3.2.1. Ring distribution. [:?] To serve as reference for the 
nongyrotropic cases to follow, in Figure 2 we show a time series for 
an unstable proton nng distribution in the perpendicular velocity 
plane. For the rest of the paper this simulation run w ill be referred 
to as the “nng case.” For the nng case we use = 10 and 
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Figure 3. Time evolution of particle diagnostics for the 
distnbution shown in Figure 2. 


= 1 .0. Figure 3 shows the corresponding time history for 
fluctuating field wave ener^^ the four dia^ostics 

defined above. Field saturation occurs at Qp/ ^ 20 in conjunction 
with df/dv^ Note the initial sharp increase in 3||p and the 

sharp decrease in (v^^) Just before saturation. While 3||p continues 
to increase, continues to decrease, but at slower rates during 
the nonlinear stage. The sharp rise in stops just after saturation. 
By this time the ring protons have completely diffused to fill in the 
center of the ring, as seen in the bottom right panel of Figure 2. The 
parameter ^>(r) starts at zero and remains at low levels throughout 
the run, as expected for a gyrotropic distribution. 

3.2.2. Nongyrotropic distribution: Run 1. [ 24 ] Section 
3.2.1 gave an overview of the effects produced by an unstable 
ring distribution. Here we begin with a similar analysis of several 
different nongyrotropic distributions. Figure 4 shows the 
perpendicular velocity space time series for a nongyrotropic 
distribution. This run uses = 1.0, = 1.0, and <bo = 

1 .0. By Op/ ^ 1 1 the particles have diffused significantly in the 
radial direction, while azimuthal diffusion is minor, and the two 
gaps are still clearly visible. Figure 5 shows the time variation of 
the four quantities defined above as well as (65/5o) for run 1 
(gray) and for the ring case (black). The field energy saturates at 
Qp/ % 17 with amplitude {hB/Bof ^ 0.04. Oscillations in the 
fluctuating field energy at jJr^Op si 0.7 are found (not shown) after 
saturation. The field energy of run I peaks sooner and at a larger 
amplitude than the ring case, while the next three parameters are 
generally the same. ^>(/) shows some correlation with the^field 
energy rise and saturation. At the first peak in {bB/Bo) the 
distribution is still strongly nongyrotropic. By the second peak in 
the field energy, the distribution has become more gyrotropic, and 


the rapid diffusion in the radial direction has slowed. An interesting 
feature to point out is the asymmetry in the radial diffusion. Figure 
4 shows that as the protons undergo their Larmor gyromotion in the 
clockwise direction about the magnetic field, the leading edge of 
the distribution in the perpendicular velocity plane diffuses faster 
than the trailing edge. This may have something to do with the fact 
that in the rest frame of the particles, the circularly polarized 
fluctuating electric and magnetic fields seem to sweep by the 
particles in a counterclockwise direction because they rotate with 
frequency ^ < Qp. Consequently, the fields encounter the leading 
edge of the distribution first. 

3.2.3. Nongyrotropic distribution: Run 2. [ 25 ] The 
distribution used for this run (not shown) is the same as in the 
previous case except the perpendicular velocity has been doubled 
to V = 2.0. The main difference between the two runs is that 
the field energy for run 2 is 6 times larger than for run 1. Recall 
that the field energy is proportional to for an unstable ring 
distribution. The results from the two runs and their comparisons 
to gyrotropic distributions are otherwise qualitatively the same. 
The main difference between these two cases is that the fields in 
run 2 show several small-amplitude spikes in the waveform while 
the ring simulation produced only a few. The spikes have 
frequencies greater than Qp. These characteristics were not seen 
in run I, implying that the purely nongyrotropic effects may only 
appear for distributions with sufficiently high perpendicular 
velocities. 

3.2.4. Nongyrotropic distribution: Run 3. [ 26 ] The next 
run is motivated by observations that indicate that nongyrotropic 
distributions could be made of multiple as well as single bunches. 
For example, Galileo observations made at a distance of -60 Re 
from Earth show an ion distribution that appears to be made of 
three separate bunches [Frank et ai, 1994]. To investigate this t^e 
of observed nongyrotropic distnbution consisting of multiple 
bunches, a simulation is run using three equidensity bunches of 
protons symmetrically spaced about the origin. All other 
parameters are the same as those used in the ring case. Recall 
that the parameter does not depend on the number of bunches, 
but only the total arc length of the removed azimuthal segment of 
the nongyrotropic velocity distribution in the perpendicular 
velocity plane. The example in run 1 therefore has the same 
value of even though the particle bunching is different than 
the three-bunch case. 

[ 27 ] The perpendicular velocity space evolution of this distri- 
bution is shown in Figure 6. The main difference between run 1 
and run 3 is that in this example, as the three bunches of particles 
gyrate around the magnetic field in the perpendicular velocity 
plane, the leading edge (head) and trailing edge (tail) of each bunch 
of particles diffuses equally in with no asymmetry. It is 
interesting that the observations of Frank et ai [1994] follow this 
same pattern; that is, the two-bunch observation shows the head- 
tail asymmetry, and the distribution with three, nearly symmetric 
bunches shows little if any asymmetry. Our simulations show 
similar results. An analytic treatment of the particle dynamics is 
required to fully understand how this asymmetry is initially 
formed. The time variation of the four diagnostics for this distri- 
bution as well as the growing electric field time histories are 
essentially the same as in run 1 . 

3.2.5. Nongyrotropic distribution: Run 4. [ 28 ] Most of the 
observations of nongyrotropic distributions cited in this text do not 
show appreciable ion background densities; however, finite 
background ion populations with gyrotropic Maxwellian 
distributions are a possibility. To investigate the effect of a 
gyrotropic background plasma on the nonlinear evolution of the 
nongyrotropic distribution, we have done several simulations that 
include a background proton population. Run 4 is one such 
simulation. Figure 7 shows the total proton distribution in the 
perpendicular velocity plane for run 4. The cold proton 
background is centered at Vy = v, = 0. For this run, v^o/v^ = 1 .0, 
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Figure 4. Velocin^ space time senes in the plane perpendicular to the background magnetic field for an initiallv cold 
nongyrotropic proton distribution with 4>o = 1.0, ' a = LO and /ing = TO. 


n„ /n^ = 0.5 (i.e.. 50®/o of the protons make up the background 
population: ^ «back = ri,). and 4>o - 1 -0 Throughout the run. 

background protons (as seen in the perpendicular velocitv' plane) are 
diffused preferentially in the direction perpendicular to the leading 
edge of the nongyrotropic proton distribution. This diffusion toward 
a Maxwell ian-1 ike distribution takes approximately twice as long as 
a comparable run with no proton background population (i.e., with 
- «c) Furthermore, (bB/Bofr^ is about half the value of the 
maximum fluctuating field energy of run 1. The effect of adding 
cold background protons to the plasma (and therefore lowering the 
number of particles in the nongyrotropic distribution) is to reduce 


the maximum fluctuatmg field energv’ in direct proportion to 
The nongyrotropic nature of the distribution remains detectable for 
a length of time inversely proportional to 

3.3. Parameter Search 

[ 29 ] The extent that fluctuating fields produced b\ nongyro- 
tropic protons differ from those fields produced b\ nng protons is 
now discussed, and the process responsible for saturation of these 
growing waves is investigated 

3.3.1. Field energy dependence on [M)] The 

dependence of the maximum fluctuating field energy on the 
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Figure 5. Time evolution of particle diagnostics for simulation 
run 1 of a nongyrotropic proton distribution, with the same 
parameters as shown in Figure 4, is shown by the gray line. The 
ring case (black line) is shown for comparison. 


nongyrotropic nature of the protons is shown most clearly in 
Figure 8. Here is plotted versus for = l.O 

(top panel) and -2.0 (bottom panel). In both examples, 

there are no background protons. As the initial proton distribution 
is made increasingly nongyrotropic (i.e., increasing ^o) the 
maximum field energy increases by up to --60% in the top panel 
and by up to ~70% in the bottom panel over that produced by a 
gyrotropic ring distribution. This result is significant because it 
clearly shows that the nongyrotropic nature of the distribution 
includes a source of free energy not found in the nng distribution. 

3J.2. Field energy dependence on perpendicular beam 
velocity. [3i] Figure 9 shows (65/fio)"max versus the perpen- 
dicular beam speed squared for a series of simulations using 
either a nongyrotropic distribution (triangles) or a ring 
distribution (circles) to model the protons. The top panel shows 
results from runs using nongyrotropic distributions with = 1 .0, 
while the bottom panel uses 4>o = 15. Observations of 
nongyrotropic distributions in the magnetotail have values of 
vjv^ within the range used in Figure 9. In all runs there are 
no background protons. The maximum field energy is seen to be 
an approximate linear function of for the runs using a 

ring distribution. For runs using a nongyrotropic distribution, this 
linear relationship holds as well, but only for low values of v^o- 
The nongyrotropic case begins to diverge from the ring case at 
large perpendicular velocities and larger 4>o» where nongyrotropic 
effects begin to dominate. 

33J. Field energy dependence on nongyrotropic proton 
density. [ 32 ] Figure 10 shows the maximum field energy 
dependence on density for ring protons (solid line) and 


nongyrotropic protons (tnangles) with = 1.5 and v_^ Va = 10.0. 
In Figure 10 a Maxwellian background proton population is used 
when ^ng ^ Observations from the magnetotail show that 
nongyrotropic distributions have low background ion densities, so 
they would fail toward the right-hand side of Figure 10. The 
nongyrotropic effects cause the maximum fluctuating energy to 
diverge sharp Iv from a linear dependence on density for 
> 0.75. 

3.3.4. Diffusion time for a nongyrotropic distribution to 
become gyrotropic. [ 33 ] Figure 11 summarizes the time it 
takes for initially nongyrotropic distributions to become gyro- 
tropic. We define the plasma to be gyrotropic when <I>(r) <0.1, 
and the time it takes for the particles to reach this point is 
denoted by (a- The top panel shows that increases with 
increasing Distributions with larger 4>„ require more 
diffusion to become gyrotropic, implying a larger /<j, but this 
longer time is partly offset because larger <l>o also imply larger 
las shown m Figure 8) and a corresponding 
increased diffusion rate of the particles. The bottom panel 
shows that decreases significantly with increasing v This 
occurs because as v is increased, fluctuating field amplitudes 
increase, causing enhanced diffusion. 

3.4. Saturation Mechanism 

[ 34 ] The time at which saturation occurs is defined here as 
the time when fluctuating field energy reaches its maximum 
value. The simulation results indicate that saturation of the 
growing fluctuating fields occurs by the same mechanism 
whether the fields are produced by a ring distnbution or by 
a nongyrotropic distribution. The saturation mechanism is 
quasi-linear diffusion of the particles leading to an increase 
in d,|p and a decrease in until these parameters approach 
the linear theory threshold values discussed by Convery and 
Gary [1997] and given by 



for ^ I • Saturation w as found to be relatively independent 
of ^o- The evidence leading to the conclusion that saturation is a 
quasi-linear process for the nongyrotropic distributions considered 
here is as follows: 

1. If the excited fluctuating field spectrum that causes the 
diffusion is sufficiently broad in frequency and/or wave 
vector, this process should be adequately described by 
quasi-linear diffusion [Coroniti^ 1985]. The linear theory 
dispersion for parallel propagating, electromagnetic, proton 
cyclotron waves in a plasma consisting of Maxwellian 
electrons and protons modeled by a cold ring distribution 
produces growing modes that are unstable over a wide 
range of wavelengths for parameters typically used in our 
simulations. Our simulation results show enough general 
similarities between waves produced by ring and nongyr- 
otropic distributions that we might expect similarities in 
their linear dispersion. 

2. In almost all of our simulations, for both ring and 
nongyrotropic distributions, saturation occurs when both 3||p 
and are close to the theoretical linear threshold condition 
given by Conxery' and Gary [1997]. This implies that quasi- 
linear diffusion, which was found to be the saturation 
mechanism for waves produced by a proton ring distribu- 
tion [Convery and Gary\ 1997], is also the cause for waves 
driven by nongyrotropic distributions. Although we have 
not carried out a linear theory study for a nongyrotropic 
plasma, we found in the simulations that nongyrotropic 
distributions saturate at about the same time as the ring 
simulations and with approximately the same values of 3|fp 








Figure 6. Velocit>' space time senes in the plane perpendicular to the background magnetic field for run 3. which 
uses an initially cold nongyrotropic proton distribution composed of three symmetncally placed bunches of particles. 
Here 4>o = 1.0, = 10, and ^ng = 1.0. 


and v\ We therefore conclude that quasi-linear diffusion 
applies in both cases 

3. Another charactenstic of quasi-linear diffusion is that the 
proton distribution should continue to diffuse toward a 
Maxwellian for some time after saturation. As shown by 
the velocity space, time history, and particle diagnostic 
plots discussed previously, the proton distribution diffuses 
rapidly in 3||p and up to saturation. The distribution then 


continues to evolve at a slightly slower rate toward a 
Maxwellian distribution at late times. 

[35] Collisions do not play a role in the saturation process 
because the Coulomb collision frequency is much lower than the 
proton gyrofrequency timescale. We rule out particle trapping 
because the proton distribution continues to diffuse toward a 
Maxwellian after saturation. If trapping played a role in saturation, 
the distribution would retain clear non >Max well lan features for 
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Figure 7. Velocity space time series in the plane perpendicular to the background magnetic field for run 4, which 
uses a plasma consisting initially of both cold nongyrotropic protons and background Maxwellian protons. Here 4>o “ 
1.0, = l.O, and ~ 0.5. Initially, the background particles are preferentially scattered along a line 

connecting the two leading edges of the nongyrotropic distribution. 


some time after saturation. Instead, in our runs, azimuthally bunched 
features continue to steadily evolve toward a gyrotropic distribution 
after saturation, with no indication of trapping during late times. 

3.5. Other Characteristics of Nongyrotropic Distributions 

[36] It has been shown that for certain parameter regimes, 
strongly nongyrotropic distnbutions can be unstable and generate 
waves that cause the distributions to quickly diffuse toward a 
gyrotropic form. This implies that it might be difficult to make 


direct observations of nongyrotropic distribution ftmctions since 
most particle instruments require at least several seconds to make a 
full three-dimensional distribution function measurement Never- 
theless, there are several observations of nongyrotropic distribution 
functions in the magnetotail and near comets as cited in this paper. 
Additionally, for the middle and distant magnetotail, there is little 
or no wave activity detected in the regions where the nongyrotropic 
distributions are observed. The question arises, if the naturally 
occurring distributions can be nongyrotropic as the observations 
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Figure 8. Maximum fluctuating field energ>' versus with (top) 
V ,,/va - I.O and (bottom) > i‘a = ^ 0. In both panels, - 

I To. Each triangle represents a separate simulation run with a 
nongyrotropic distribution. The dotted lines connecting the 
tnangies are used simply to clanfy the trend of increasing 
maximum field energv with increasing 


suggest, w’h> are the> not unstable and diffuse toward a more 
gyrotropic distnbution'.' We address these important questions in 
this section and show that there are certain parameter regimes in 
which nongyrotropic distributions are stable Specifically, we 
examine the effect of including a finite perpendicular thermal 
spread to the ring and nongyrotropic distributions. We also inves- 
tigate single bunched distributions as opposed to the pre\ ious 
sections, where we considered distributions made of two or three 
symmetrically spaced bunches. 

3.5.1. Finite perpendicular thermal spread. [ 37 ] Observed 
nongyrotropic distributions generally show' a finite thermal spread in 
both the parallel and perpendicular velocities, that is. both Ty ^ 0 and 
T ^ 0. We have, to this point, used model ring and nongyrotropic 
distributions that ha\ e a finite parallel temperature but are cold in the 
perpendicular direction, that is, = 0. We now include a 
perpendicular thermal spread in our model nongyrotropic 
distribution. To do this, the perpendicular velocity (v^^) of the 
nng protons is modified by adding a perpendicular thermal speed 
vth ^ such that r p = Tijp. All examples discussed in this subsection 
wilFuse Tjip = T^p. For relatively low perpendicular thermal spreads, 
the simulation results are very' similar to the cold case. As the initial 


temperature is increased, however, the nongyrotropic distribution 
eventually becomes stable, independent of the value of 4>o 

[38] The perpendicular velocity space time series for this non- 
gyrotropic case is shown in Figure 12, and the particle diagnostics 
are given in Figure 13. The effects of using a perpendicular thermal 
spread in either the initial ring distribution or nongyrotropic 
distribution, relative to runs using cold initial distributions, are 
twofold. First, it reduces the growth rate of the instability' so that 
the non-Maxwellian features are more slowly diffused away, that 
is, rip/ ~ 40 for this example compared to ^ 25 for the 
otherwise identical cold case. Second, the maximum fluctuating 
field energy is slightly less than in the cold case. Overall, however, 
the results are qualitatively similar to the cold case. 

[39] One similarity that the cold and thermal runs using non- 
gyrotropic distributions have is that in making the distribution 
more nongyrotropic, (65 -5o)max increased over the correspond- 
ing gyrotropic ring case. Tfiese results provide an answ'er to the 
question about the effects of allowing a finite perpendicular 
temperature in the nongyrotropic distribution. Namely, w hile w arm 
distributions are a closer representation of the observed distribu- 
tions, the similarity of our cold and warm results validates our 
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Figure 9. Maximum fluctuating field energy versus perpendi- 
cular beam speed squared for a series of simulations using (top) 4>c, 
= 1.0 and (bottom) <I>o = 15. In both panels, n^^ric = 10. Each 
p>oint represents a separate simulation run. The circles indicate that 
a ring distribution was used in the simulation, and a triangle 
denotes that a nongyrotropic distribution was used. All other 
parameters are the same 
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Figure 10. Maximum fluctuating field energy versus ^ng for 
simulations using <Do = 1.5 and v = 10.0. Each triangle 
represents a separate simulation using a nongyrotropic proton 
distribution. The solid line is the maximum fluctuating field energy 
due to an equivalent, cold gyrotropic ring distnbution. 


[42] Two important conclusions can be drawn from these 
results. First, model distributions with a finite thermal spread in 
give results qualitatively similar to simulations using initially 
cold perpendicular velocity. However, there are small quantitative 
differences between the results depending on whether the cold or 
thermal distnbution is used. Therefore the thermal model, which is 
a more accurate representation of observed nongyrotropic magne- 
totail plasmas, should be used to make detailed quantitative 
compansons with magnetotail observations of quantities such as 
fluctuating field amplitude and saturation times. Second, for the 
particular model and limited parameter regime used here, non- 
gyrotropicity (i.e., ^ 0) is not a sufficient condition for the 

production of fluctuating fields: that is, making a stable nng 
distnbution nongyrotropic does not necessanly produce an unsta- 
ble plasma. These results indicate that the nongyrotropic nature of 
a plasma population alone does not generate significant fluctuating 
fields and therefore may be considered stable for some parameters 
regimes even though the nongyrotropic distribution is non-Max- 
wellian. 

3,5.2. Single bunch distributions. [ 43 ] Up to this point our 
model nongyrotropic distributions have all had one property in 
common. In the perpendicular velocity plane the nongyrotropic 
particles are confined in gyrophase in groups of equal number 
density that are spaced symmetrically about the origin. This model 
IS used because the structure of most hybrid codes, including ours. 


initial conclusions gathered from simulations using an initially cold 
{T_p = 0) model distnbution. 

[ 40 J Is a distnbutiqp function with a finite nongyrotropic nature 
(^o ^ 0) a sufficient condition to produce instabilities, or does it 
simply enhance the field energy of a gyrotropic plasma which is 
already extant in an unstable parameter regime? This is the second 
question we seek to answer in this section, that is, is there any 
parameter regime in which a gyrotropic ring distribution is stable 
and yet an otherwise identical distribution, but with a nongyro- 
tropic nature (^o ^ A), produces an instability? 

[4t] Our studies show that increasing the initial ring perpendic- 
ular thermal spread further produces a decrease in the maximum 
fluctuating field energy. At a temperature corresponding to 3i|p = 
0.5 (with r^p = 71|p), the ring distribution is stable. To examine this, 
a run was carried out using a nongyrotropic distribution created by 
removing two 45® sections (4>o = 15) from the thermal ring 
distribution (not shown). The protons in this case did not diffuse 
significantly until Qpt % 50. The maximum fluctuating field energy 
reached a peak wave energy that was only ~1% of the energy 
produced by an identical run using a cold initial distribution. There 
was little activity in the other diagnostics except for the slow, 
gradual decrease in 4>(/) with time. Additionally, two other runs 
(not shown) were done using warm nongyrotropic distributions 
with parameters identical to the first warm nongyrotropic distribu- 
tion case except for changing to 0.5 and 1 .0, respectively. In all 
three cases the nongyrotropic distribution remained relatively intact 
with low wave activity, independent of the value of 4><j. These 
examples imply that 0 is not a general condition for 

instability for this parameter regime. It must be linked to the 
requirement that the nongyrotropic plasma would already be 
unstable in a related gyrotropic manifestation, at least for thermal 
parameter values studied here. For unstable distributions, however, 
works well as a diagnostic to encompass the effect of a 
distribution s nongyrotropic nature on the magnitude of the max- 
imum fluctuating field energy. It is important to point out that in 
certain parameter regimes other that what is considered here, the 
introduction of a nongyrotropic plasma distribution by itself can 
generate an instability when the corresponding gyrotropic situation 
is stable. This is true for both stationary [Brinca and Romeiras, 
1998; Brinca et aL, 1998; Brinca. 1999; Motschmann and Glass- 
meier. 1998; Motschmann et al., 1999] and time-varying [Motsch- 
mann et ai. 1999] nongyrotropic distributions. 
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Figure 11. This shows the time for the initially nongyrotropic 
distribution to become gyrotropic (to) versus (top) <l>o and (bottom) 
v^o, where Iq is defined as the time it takes for the plasma to reach 
^(0 < 0 . 1 . 
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Figure 12. Velocm space time senes in the plane perpendicular to the background magnetic field for a 
nongyrotropic proton distribution with an initial thermal distribution of particles about v Here “ 1 ' ' a = 

1 .0, ~ 1 -O' ~ 0 1 


requires us to assume that there are no net. zeroth order currents. 
While man\ insightful results have been achieved by using this 
symmetric, nongyrotropic model, it is worthwhile to explore 
nonsymmetric distributions because the observed nongyrotropic 
distributions are typically not symmetrically spaced about the 
magnetic field. Second, the question arises as to whether the 
conclusions reached ui the analysis of symmetric nongyrotropic 
distributions apply to nonsymmetric, nongyrotropic distributions as 
well. Up to now' we have parameterized all nongyrotropic effects in 
a single term, that does not take into account the number or the 


spacing of bunches in a given distribution, only their total a 2 umuthal 
spread in Ao. For example. 4^^ has the same value for a distribution 
made up of two 90*^ bunches or a single 1 80^ bunch Is 4^^ really a 
valid way to parameterize the nongyrotropic nature of a 
distribution'!^ These issues are examined in this section. 

[ 44 ] Two asymmetnc nongyrotropic distributions with differ- 
ent values of 4>e are studied in this section. Relaxing the 
requirement of symmetry in perpendicular velocity' requires the 
inclusion of a finite proton background density with an offset 
velocity set to maintain a net zero initial current Here we are 
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Figure 13. Time evolution of particle diagnostics for the 
nongyrotropic distribution shown in Figure 12. 


concerned mainly with the qualitative aspects of single bunched 
distributions. In general, larger nongyrotropic densities would 
most likely produce larger fluctuating field energies for a given 
value of v^q/v^. We have been using = l.O throughout 

much of this paper; however, we know that a finite background 
density lowers the maximum field energy (see Figiue 10), and 
we have found that the simulations of this section do not show 
appreciable growth when using = 1.0. At ^ 2.0 

we begin to see some wave growth; that is, {bB/Bo)max ~ 0 01. 
For all the nongyrotropic examples considered in this section we 
use vxo/va = 3.0 and = 0.2, and for comparison, a run 

using a cold ring with = 0.2 and = 3.0 produced 

growing field energy with a peak value of (65/5o)max ^ 0.03. 
At saturation (Qp/ ^ 20) the ring distribution was nearly 
Maxwellian owing to its diffusion in both the parallel and 
perpendicular velocity. 

3.5.2. 1. Single bunch nongyrotropic distribution: Example 

1; [45] This run, shown in Figure 14, uses a nongyrotropic 
distribution created by removing a 90° section ~ 0.5) from a 
cold ring. The field energy reaches a peak value of (65/.5o)^ ~ 
0.035 at Qpt =55 16. Figure 14 shows that the distribution still retains 
some signs of its initial nongyrotropic character at this time. The 
main result here is that this single bunch distribution is unstable 
and produces fluctuating field energy larger than in the ring case. 

3.5.2.2. Single bunch nongyrotropic distribution: Example 
2: [ 46 ] In this example a 270° section (4>o = 1 .5) is removed from 
a cold ring. Figure 15 shows the perpendicular velocity space time 
history of this distribution. The field energy reaches a peak value of 
(&5/5o)^ ~ 0.06 at Qpi ^ 18. Thus the maximum field energy, 
compared to the ring case and to the previous example I , increases 
with increasing 4>o- 


[47] From the above examples we conclude that single bunch 
distributions are similar to the nongyrotropic distributions made of 
two or more symmetncally spaced bunches that were examined in 
the earlier sections of this paper. Their characteristics include an 
increase in maximum field energy with increasing and 
although significant diffusion has taken place in v: and by the 
time of saturation, the distributions are still nongyrotropic. The 
main result of this subsection is that the stability of nongyrotropic 
distributions does not seem to depend on the detailed azimuthal 
spacing or number of bunches that make up the distnbution. This is 
relevant because nongyrotropic distributions with any number of 
bunches have been observed in the Earth's magnetotail. We can 
now have more confidence that the parameter which does not 
take into account the number of gyrobunches. gives an accurate 
descnption of the maximum field energy increase due to non- 
gyrotropic effects. 

[ 48 ] The goal of this section has been to determine which 

parameters, if any, stabilize nongyrotropic distnbutions. To this 
end, we looked at distributions with finite perpendicular thermal 
spread and those with single bunched distnbutions. First, we 
isolated the nongyrotropic nature of the protons by using a 
distnbution characterized only by 4>o — 0 (i.e., <7/ > 0) so 

there was no temperature anisotropy or any other free energy 
sources, such as nng-like {dfidv _ > 0) effects, associated with the 
plasma. We found that large values of r_p can cause such a 
nongyrotropic distribution to be stable, even for very large values 
of Iti this way it is shown that 7^ 0 is not necessarily a 
sufficient condition for instability in some parameter regimes. It 
is more accurate to say that if a gyrotropic distnbution, such as a 
ring, is unstable, then nongyrotropic modifications to the plasma 
(as quantified by 4>o) will enhance the growing fluctuating field 
energy. In such cases, accurately descnbes the enhanced 
fluctuating field energy due to the nongyrotropic nature of the 
plasma. A new parameter that can describe the effect of the 
nongyrotropic nature of the plasma in all parameter regimes 
should be developed in the future. This new parameter should 
be a function of and another variable that accounts for the 
gyrotropic stability of the plasma. 

[49] We summarize that nongyrotropic plasmas may be stable if 
they are characterized by high temperatures in the parallel and/or 
perpendicular velocity. Other factors contributing to stability are 
the proton density ratio and the perpendicular beam speed; that is, 
stability is likely if the nongyrotropic (or ring) distribution density 
is a sufficiently small fraction of the total plasma density or the 
particles have small enough values of v o.'va such that linear 
stability threshold of an otherwise identical gyrotropic ring dis- 
tribution [Conver}' and Gary^ 1997] cannot be crossed. 


4. Conclusions 

[50] This paper has explored the stability and behavior of 
nongyrotropic proton distributions. It was motivated by observa- 
tions from the Geotail [Saito et al., 1994] and Galileo [Frank et ai, 
1994] spacecraft revealing nongyrotropic distnbutions in the 
Earth’s magnetotail. These data raised the possibility that standard 
gyrotropic, Maxwellian-like modeling of particle distributions is an 
inappropriate description of the plasma in some regions of the 
Earth's magnetosphere. 

[5!] Our nongyrotropic model distribution has allowed us to 
determine consequences of nongyrotropic space plasmas through a 
comparison with the effects of ring distributions. We conclude that 
the maximum electromagnetic fluctuating field energy increases as 
the initial proton distribution is made more nongyrotropic. For the 
parameters used here, this increase can be as much as 2 times the 
energy produced by the case using an otherwise identical gyro- 
tropic distribution. The maximum electromagnetic fluctuating field 
energy for a given ^^so increases with increasing and 

with increasing The fields stabilize by diffusion in the 
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Figure 14. V^locit> space time senes in the plane perpendicular to the background magnetic field for the case with 
an asymmetnc. non^otropic distribution function w'lth a three-quarters ring. 


parallel and perpendicular velocity, leading to an increase m 3|jp 
and a decrease in to values approaching the linear theor>' 
instability threshold of Convery and Gan [1997]. The saturation 
mechanisms and the timing are relatively independent of 
time for the nongyrotropic distribution to diffuse in phase angle 
toward a g>TOtropic distribution {to) is found to increase with 4>o 
and decrease with v^^/va- This agrees with Cao et al. [2000]. who 
showed that pickup ions in solar wind plasmas have increased 
phase angle diffusion with increased injection velocity'. 

[52] Cao et al. [1995] and Motschmann and Glassmeier [1993] 
studied the nongyrotropic nng instability near a cometaiy environ* 


mcnt and found it to be linearly unstable. The linear growth rate of 
this instability was found to be greatest for a ring distribution and 
to decrease with an increasing nongyrotropic nature of the distri- 
bution. The simulation results described in this study show that 
ibB increases with the nongyrotropic naUire (4>o) of the 

distribution function, as demonstrated in Figure 8. In comparing 
these two results, it should be piointed out that the types of 
nongyrotropic distributions considered by Cao et al. [1995] have 
a different form than the model used here and the parameter 
regimes of both Cao et al. [1995] and Modchmann and Glass- 
meier [1993] differ from ours. Thus the linear gro^^th rates cannot 
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[ 54 ] When the gyrotropic plasma is stable, as in the case for a 
sufficiently warm ring, making the gyrotropic plasma nongyro- 
tropic does not produce instability for the parameters used in this 
study. When the plasma is in an unstable parameter regime, 
nongyrotropic distributions produce growing electromagnetic pro- 
ton cyclotron w aves similar to those produced by gyrotropic proton 
ring distributions, and increasing the nongyrotropic nature (i.e., 
increasing ^o) will enhance the maximum fluctuating field energ>. 
The growing modes in both the ring and the nongyrotropic cases 
are parallel propagating and occur at low ^ Op) frequencies 
These results lead to the conclusion that the introduction of a 
nongyrotropic modification to the ring distribution enhances the 
fluctuating fields that would be produced by a ring distribution. 
The nongyrotropic nature of a plasma population does not by itself 
necessarily ensure that the plasma is unstable, at least for the 
parameters studied here. 

[55] Similar conclusions were made on the basis of linear 
stability studies by Brinca and Romeiras [1998] and Brinca 
[1999], who showed that in certain parameter regimes, parallel 
propagating waves are stable in the presence of time-varying 
nongyrotropic distributions, unless a finite perpendicular current 
exists. Moischmann et ai. [1999] found that a time-varying non- 
gyrotropic plasma is linearly stable to parallel propagating waves, 
for certain parameters, unless an additional energy contribution 
such as a temperature anisotropy is present. Motschmann et ai 
[1999] also states that nongyrotropic distributions are expected to 
be stable when both the perpendicular and parallel velocities are 
conserved and particle diffusion occurs along the gyrophase angle 
only. This is because particle energy is independent of phase angle 
When the gyrophase diffusion is in conjunction with perpendicular 
and/or parallel velocity' diffusion, the particle energy can be 
converted to wave energy, resulting in growing waves [Motsch- 
mann et ai, 1999]. 

[56] A difference found in some of the electric field time 
histories is that the nongyrotropic distribution is associated with 
waves having intermediate, fast spikes in amplitude in addition to 
the low -frequency component near ~ flp. The low-frequency 
component is similar to that seen in the simulations using a ring 
distribution. The nongyrotropic distribution produces these waves 
w^ith little or no linear growth phase, while the ring case shows 
clear linear growth. 

[57] We have discussed observational evidence for the existence 
of nongyrotropic distributions in the Earth’s magnetotail and have 
established that these distributions lead to enhanced parallel 
propagating electromagnetic waves. The question that anses now 
is, what are the magneiospheric consequences of these fluctuating 
fields? Our simulations show that the growing electromagnetic 
fields produced bv the nng and nongyrotropic distributions saturate 
by difhision of the particles in both the parallel and perpendicular 
directions. The considerable diffusion in the parallel direction 
could contribute to particle entry' into the loss cone, thereby leading 
to enhanced precipitation. This is significant since particle precip- 
itation calculations that employ the concept of a loss cone tradi- 
tionally use distributions that are azimuthallv symmetric. If 
fluctuating fields produced by azimuthallv bunched distributions 
change the rate at which particles are scattered into the loss cone 
compared to a gyrotropic distribution, then it will be necessary for 
these calculations to be revised for some situations to include 
nongyrotropic effects. Our simulations indicate that nongyrotropic 
particles are associated with fluctuating fields having enhanced 
wave energy above that of the gyrotropic case. This leads to an 
increase in the diffusion rate over the gyrotropic case 

[58] If field lines in the outer boundary' of the magnelotail 
plasma sheet are connected to the auroral oval, this indicates that 
particles lost in the w'ay described above could contribute to the 
aurora and low-altitude ion flux coming from the tail. This loss 
mechanism is significant because it affects particles in the central 
plasma sheet that otherwise would be unlikelv to have V| large 


enough to enter the loss cone. For example, large-scale kinetic 
particle trajectory studies show' that in the absence of w av e-particle 
interactions, only a small fraction of plasma sheet boundary layer 
ions have pitch angles that are small enough to be able to 
precipitate into the auroral ionosphere [Ashour-Ahdalla et ai, 
1993a]. There have been other suggestions for maintaining ion 
precipitation from this region. For example, Zelenvi et ai [ 1 990] 
suggested that stochasticity' could lead to an effective diffusion in 
the CPS resulting from particles undergoing many crossings of the 
field reversal region and their subsequent acceleration. Neither this 
mechanism nor the process discussed here of increased ion 
precipitation due to instabilities associated with nongyrotropic 
distributions rules out the other. These processes both involve 
acceleration in the magnetotail field reversal region, and thev may 
be related or may be operating simultaneously to effect the 
observed ion precipitation. The results here indicate that the effects 
of waves in the presence of nongyrotropic distribution functions 
must be considered self-con si stently to fully understand ion pre- 
cipitation from the magnetotail. 
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